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We study the physical properties of ZnX (A=0, S, Se, Te) and CdX (A=0, S, Se, Te) in the zinc- 
blende, rock-salt, and wurtzite structures using the recently developed fully ab initio pseudo-hybrid 
Hubbard density functional ACBNO. We find that both the electronic and vibrational properties of 
these wide-band gap semiconductors are systematically improved over the PBE values and reproduce 
closely the experimental measurements. Similar accuracy is found for the structural parameters, 
especially the bulk modulus. ACBNO results compare well with hybrid functional calculations at a 
fraction of the computational cost. 


I. INTRODUCTION 

Among II-VI semiconductors, oxides and chalcogenides 
with Cd and Zn are the subject of extensive scrutiny 
for their potential applications in spintronics, optoelec¬ 
tronics, and photovoltaics.lll These applications rely on 
hetero-junctions, quantum dots, and other nanostructures 
where the delicate interplay between structural and elec¬ 
tronic properties must be carefully captured to understand 
and improve the device’s performance. Common compu¬ 
tational approaches based on density functional theory 
(DFT), however, fall short when trying to predict accu¬ 
rate structural and electronic properties, greatly limiting 
the development process. 

In DFT, most of the calculations involving II-VI semi¬ 
conductors and oxides use the exchange-correlation (xc) 
potential within the local density approximation (LDA) or 
the generalized gradient approximation (GGA).l^ While 
the DFT-LDA/GGA method provides reasonable predic¬ 
tions for structural properties, it fails dramatically in de¬ 
scribing the electronic properties, especially the band- 
gaps, which are underestimated by 50% or more with re¬ 
spect to experimental values. In ZnO, for instance, the cal¬ 
culated band gap is 0.8 eV, while the experimental value 
is 3.44 eVj^GdO is described as a semi-met al,l^ whereas 
experimentally rocksalt GdO has an indirect band gap of 
0.9 eV.I^ Similar discrepancies in the calculated and mea¬ 
sured band-gaps are found in the chalcogenides as well.l^ 
In the case of Gd- and Zn- oxides and chalcogenides, in 
addition, the energy of the occupied d-manifold is found 
about 3 eV higher than the experimental value and lead 
to fictitious d-p mixing in the valence band.^^ 

The limitations mentioned above are a consequence of 
the lack of self-interaction corrections^ and of the deriva¬ 
tive discontinuity in the exchange-correlation energ}®®^ 
in all the traditional functionals such as LDA or PBE- 
GGA. Glearly, only the inclusion of non-local correlations 
and non-local Hartree-Fock (HF) exchange can, in part, al¬ 
leviate these problems. A number of approaches have been 


proposed to overcome these deficiencies, among which: 
LDA plus self-interaction (LDA-SIC),^ self- intera ction- 
relaxation correction (SIRC-LDA),1^ DFT+Uj ^^ * ^^ * differ¬ 
ent version s of the GW approximation,!!^®^ and hybrid 
functionals.!!^!^ Among these approaches, the ones that 
h ave g ained more traction in recen t years are the DFT+ 
LP!^ and hybrid functionals .!!^^^ 

In DFT+U one introduces on-site Goulomb and ex¬ 
change interaction terms, U and J, to account for the 
localization of d states. The on-site orbital-dependent pa¬ 
rameters U and J correspond to the Coulomb and ex¬ 
change couplings between electrons of a particular angu¬ 
lar momentum that are localized on the same atom. The 
appeal of DFT+U relies on its effectiveness and low com¬ 
putational cost in correcting for the over-delocalization of 
the d electrons in transition metal ions. Although the orig¬ 
inal DFT+U formulation is rigorous, in most cases U and 
J are treated as empirical adjustable parameters which are 
often obtained by fitting the band structure to available 
experimental values.!®! This approach requires fitting in¬ 
formation from experimental data and its predictive value 
for new materials and hetero-structures is limited. There 
are a few common ah-initio methods to derive U and J 
such as t he con strained LDA (cLDA)[®!and linear response 
approachThese methods are, unfortunately, compu¬ 
tationally expensive, often requiring large supercell calcu¬ 
lations. Moreover, although the linear response method 
has been widely used for open-shell systems, it is not suit¬ 
able for closed-shell elements such as Zn and Gd, where 
the localized bands are completely full and insensitive to 
linear perturbations.!®! 

Hybrid functionals are based on the idea of computing 
the exact exchange energy from the Kohn-Sham wavefunc- 
tions and to mix it with the (semi)local approximation of 
correlation energy of DFT.!®! this respect, the method 
does not suffer of any limitations in dealing with closed- 
shell elements and it is very successful in predicting the 
energy gap for semiconductors and insulators. Clearly, 
some degree of exact exchange is necessary for a more ac- 
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curate description of the electronic structure. However, 
also hybrid functionals have a somewhat empirical com¬ 
ponent, since the level of mixing is not determined from 
first principles.!^ 

In order to facilitate the accurate characterization of 
electronic properties of materials at a low computational 
cost, a fundamental condition for the development of ef¬ 
fective high-throughput quantum-mechanics frameworks 
for accelerated materials discovery,!^^^^^ some of the au¬ 
thors have recently introduced the ACBNO functional,!^ 
a pseudo-hybrid Hubbard density functional that intro¬ 
duces a new ab-initio approach to compute U and J that 
does not contain any empirical parameter. In this work we 
demonstrate that, by using ACBNO, we can improve sub¬ 
stantially the predictive value of DFT calculations for the 
physical properties of Zn and Cd oxides and chalcogenides 
at a fraction of the computational cost of hybrid func¬ 
tionals. The paper is organized as follows: in Sec. [n| we 
briefly discuss the ACBNO method. In Sec. Ell we discuss 
our results for the Zn and Cd chalcogenides by compar¬ 
ing structural, electronic, and vibrational properties with 
more standard approaches. 

II. METHODOLOGY 

ACBNO is based on the DFT+t/ energy functional as 
formulated by Dudarev (Ref. l3T]) : 

EbFT-\-U = ^DFT + Eu — E]jc 

where F^dft is the DFT energy calculated using a LDA 
or GGA functional, the two parameters U and J have 
been replaced by an effective on-site Coulomb interaction 
^eff = U — J and E]jc takes care of the double counting 
terms in the energy expansion. For technical details on 
the ACBNO formalism we refer the reader to Ref. |30l 

ACBNO resolves the ambiguities in DFT+U by comput¬ 
ing on-the-fly the local Coulomb {U) and exchange (J) 
integrals for the specific orbitals under consideration via 
a procedure based on an ad-hoc renormalization of the 
density matrix. In this way, the value of I/eff results a 
functional of the electron density and depends directly on 
the chemical environment and the crystalline field. 

In traditional formulation of DFT+U,!^the correction 
term U was based explicitly on the localization of the d 
orbitals and was not used for the p or 5 orbitals which 
tend to be less localized compared to d electrons. How¬ 
ever, in the ACBNO formulation I/eff is computed directly 
on the chosen Hubbard center from the Coulomb and ex¬ 
change Hartree-Fock energies of the solid and can be eval¬ 
uated for any contributing orbital. The evaluation of I/eff 
for any orbital symmetry is not a new concept (see for 
instance Ref. [32j) but it is particularly relevant for the sys¬ 
tems studied in this work. 

In the current implementation of ACBNO, I/eff is eval¬ 
uated through a self consistent procedure where we start 
by calculating the electronic structure for an initial guess 
of = 0 eV and converge to I/eff within 10“^ 


eV. A plot of the convergence of I/eff with the iteration 
steps is provided in the Supplemental Information. Tab. 
[T|and Tab. [TT] provide the converged I/eff values for all the 
semiconductors studied in this work. 

ACBNO calculations have been done using the 
Quantum EspresscP^ and WanT!^^!^ packages with 
norm-conserving pseudopotentials from the pslibraryl.O 
database. A kinetic energy cut-off of 350 Ry and a k-point 
mesh of 12x12x12 was used for all total-energy calcula¬ 
tions. Hybrid functional DFT calculations were performed 
using the Vienna ab-initio simulation package (VASPj^ 
with the HSE06.!^This functional is defined by replacing 
25% of the PBE exchange interaction by a screened non¬ 
local functional with an inverse screening length of 0.2/A. 
A 6x6x6 Monkhorst-Pack k-point mesh was used and a 
plane-wave cut-off of 500 eV was used for all the HSE06 
calculations. 


III. RESULTS 

This section discusses the structural, electronic, and vi¬ 
brational properties of ZnO, ZnS, ZnSe, ZnTe, CdO, CdS, 
CdSe, CdTe computed with the ACBNO functional. 

Under standard conditions ZnO is stable in the wurtzite 
(wz) structure,!!! CdO in the rocksalt (rs) structure, while 
the rest of the compounds are stable in the zinc-blende 
{zb) structure. Due to the importance of these materials 
in semiconductor nano-structures such as heterojunctions 
and quantum dots that are usually grown epitaxially with 
far-from-the-equilibrium techniques, we study all the three 
competing phases for all the chemical compositions. These 
different structures can be stabilized by alloying and epi¬ 
taxial strain. Results not reported in the manuscript can 
be found in the Supplemental Information. 


A. Structural properties 

The lattice constants (ao) and bulk moduli (B) for each 
of the H-VI semiconductors in the three bulk phases are 
obtained by fitting the total energy as a function of volume 
to the Murnaghan’s equation of state (eos). As discussed 
in Sec.jHj Ueff is a functional of the electron density and is 
thus dependent on the geometry of the ground state. For 
a small volume change around equilibrium, the calculated 
Ueff varies linearly with the volume. As an example. In 
the top of Fig. we have plotted the variation of Ueff for 
CdO in the rs phase where we observe the linear scaling 
for both the Cd’s and the O’s Ueff, most likely related to 
bond length changes that affect the wave-functions overlap 
and the electron localization. 

In Tables. |TTI| and [IVl we have listed all the values of the 
computed lattice constants and bulk moduli of all the H-VI 
semiconductors in the three phases using PBE, HSE and 
ACBNO (and SIRC, whenever available) with references to 
existing experimental values. Experimentally, the lattice 
constants follow the trend of ZnO < ZnS <ZnSe < ZnTe 
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TABLE L Converged values of Leff (in eV) for the Zn 3d and 
the anion p states. The calculations were performed at the 
equilibrium lattice constants in Table [Till 


Phase 

ZnSd 

Anion p 

ZnO wz 

13.19 

5.57 

zb 

13.24 

6.02 

rs 

13.52 

6.05 

ZnS wz 

13.48 

3.50 

zb 

13.21 

3.42 

rs 

14.54 

3.76 

ZnSe wz 

15.00 

2.55 

zb 

14.60 

2.52 

rs 

15.68 

2.77 

ZnTe wz 

16.00 

2.40 

zb 

16.89 

2.43 

rs 

17.08 

2.66 


TABLE IL Converged values of t/eff (in eV) for the Cd 4d and 
the anion p states. The values are for the equilibrium lattice 
constants in Table imi 


Phase 

Cd 4d 

Anion p 

CdO wz 

10.69 

4.50 

zb 

10.49 

4.03 

rs 

10.73 

3.92 

CdS wz 

11.07 

3.69 

zb 

10.89 

3.64 

rs 

11.77 

3.94 

CdSe wz 

11.58 

2.68 

zb 

11.69 

2.70 

rs 

12.51 

3.05 

CdTe wz 

12.95 

2.57 

zb 

13.00 

2.55 

rs 

13.70 

2.46 


and similarly for the CdX series. This trend of increasing 
lattice constant is described in the PBE, HSE and ACBNO 
functionals. The experimental bulk moduli shows a reverse 
trend Le., the bulk modulus decreases as we go from O to 
S to Se to Te. 

In Figures [^a) and[^b), we have plotted the relative 
percentage error of the three functionals PBE, HSE, and 
ACBNO in predicting the lattice constants and bulk mod¬ 
uli with respect to the measured values. In the plot, we 
have included only the stable phases for each of the eight 
semiconductors since we did not have experimental data 
on the non-equilibrium phases. The values for the other 
three phases for the II-VI semiconductors are listed in Ta¬ 
ble |in[ Except for CdO and ZnO, all the other semicon¬ 
ductors have a stable zb structure. For ZnO, we have 
listed the error for one of the lattice parameters, a, while 
for CdO we have used the rs phase. 

The lattice constants and bulk modulii are more accu¬ 
rately described by HSE and ACBNO compared to PBE 
functional calculations. The PBE overestimates the lattice 
constants significantly (up to 2%) as seen in Fig. [^a) and 
severely underestimates the bulk-moduli (See Fig. [^b)); 
in some cases the error is more than 50% compared to 
the experimental value. The HSE functional shows a bet¬ 
ter improvement in the prediction of lattice constants re¬ 
ducing the error to less than 1%. The lattice constants 
predicted by the ACBNO follows the HSE closely with an 
error less than 1%. Note that the ACBNO preserves the 
accuracy of the HSE calculations with 8-10 times speed 
up in terms of calculation time. 



FIG. 1. (Color online) Lower panel: Energy-volume curve of 
CdO in the rs phase as calculated by the ACBNO functional. 
Upper panel: Variation of the converged Ueff of Cd and O in 
CdO as a function of lattice constant. 


The bulk modulus predicted by ACBNO shows a tremen¬ 
dous improvement over both the HSE and the PBE func¬ 
tional. In CdS, CdSe and CdTe, the agreement with ex¬ 
perimental values is exceptional with an error less than 
(0.5%). Note that the PBE error in the bulk moduli 
are greater than 20%. These results indicates that using 
ACBNO for the calculations of phonon spectra or any lat¬ 
tice dynamical property will yield much improved results, 
as already observed in the case of ZnO,where the values 
of the Zn and O Ueff’s, originally fitted to reproduce the 
experimental band gap and position of the d bands, are in 
remarkable agreement with the ones predicted by ACBNO. 
Incidentally, Ueff values fitted to HSE calcul ation s have 
been proposed in the past for ZnO and ZnS.I^^^^ While 
for ZnO the values are again close to our predictions, the 
values for ZnS are very much different. See Sec. |HI C| for 
a comprehensive discussion of the lattice dynamical prop¬ 
erties of these systems. 


For all other phases, the ACBNO functional reproduces 
the HSE values closely as seen in Table and |IV| at a 
much lower computational cost. Overall the the structural 
properties of the H-VI semiconductors computed with the 
ACBNO functional show better agreement with respect to 
experimental results. 
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TABLE III. Comparison of the lattice constants for the different compounds calculated using ACBNO and other functionals. All 
the PBE and HSE values are computed in this work. Eor wz^ the first value is the in-plane lattice parameter a and the second 
value is the out-of-plane lattice parameter c. All the SIRC values reported here are taken from Ref. [8] 


System 

Phase 

PBE 

HSE 

SIRC 

ACBNO 

Experiment 

ZnO 

wz 

3.283, 5.309 

3.260, 5.221 

3.290, 5.29 

3.270, 5.164 

3.258, 5.22(P 


zb 

4.670 

4.582 


4.580 

4 62()41J421 


rs 

4.370 

4.278 


4.289 


ZnS 

wz 

3.880,6.300 

3.850,6.271 

3.830,6.280 

3.851, 6.278 

3.811, 6.23#^ 


zb 

5.489 

5.432 

5.421 

5.437 

5.41(P 


rs 

5.143 

5.077 


5.080 

5 06()Ml45l 

ZnSe 

wz 

4.043,6.70 

4.030,6.620 


4.020, 6.630 

3.980, 6.53(F^ 


zb 

5.771 

5.708 


5.696 

5 gg7 1401461 


rs 

5.401 

5.330 


5.321 


ZnTe 

wz 

4.366,7.176 

4.35,7.132 


4.340, 7.140 

4.320, 7.10(P 


zb 

6.199 

6.154 


6.148 

g 08#I3S1 


rs 

5.788 

5.743 


5.732 


CdO 

wz 

3.700,5.850 

3.650,5.800 


3.590,5.750 



zb 

5.172 

5.103 


5.054 



rs 

4.771 

4.722 

4.781 

4.689 

4 gQ^6|48|49| 

CdS 

wz 

4.206,6.85 

4.223,6.951 

4.154,6.762 

4.183,6.682 

4.135,6.701'^® 


zb 

5.963 

5.901 

5.851 

5.892 

5 82(J40|44|51| 


rs 

5.890 

5.466 


5.459 


CdSe 

wz 

4.39,7.1 

4.350,7.112 

4.292,7.021 

4.33,7.072 

4.310,7.01(P 


zb 

6.239 

6.158 

6.071 

6.142 

6.08#' 


rs 

5.800 

5.697 


5.680 


CdTe 

wz 

4.550, 7.451 

4.580,7.460 


4.600, 7.498 



zb 

6.621 

6.571 

6.401 

6.560 

6.48(P' 


rs 

6.130 

6.092 


6.076 



TABLE IV. Comparison of the bulk moduli B (GPa) of different compounds using three different functionals. All the PBE and 
HSE values are calculated in this work. 


System 

Phase 

PBE 

HSE 

SIRC 

ACBNO 

Experiment 

ZnO 

wz 

127 

143 

159 

157 

136-183»i4ui4a 


zb 

123 

143 


162 



rs 

156 

188 


205 

177-22^^ 

ZnS 

wz 

60 

74 

81 

79 

7#^ 


zb 

68 

74.4 

81 

80 

76.#^ 


rs 

84 

95.4 


101.6 


ZnSe 

wz 

57 

62.69 


64 



zb 

55.8 

61 


66.3 


rs 

70 

78 


66.8 


ZnTe 

wz 

45 

58.46 


56 



zb 

43 

47.19 


49.4 

50.d® 


rs 

53.3 

58.82 


60 


CdO 

wz 

92.7 

104 


114 



zb 

91 

102 


124 



rs 

119 

137 

152 

170 

14’;E11 

CdS 

wz 

54 

57 

74 

53 

610 6^ 


zb 

53.5 

58.94 

70 

63.3 

5P1 


rs 

63.3 

78.6 


83.2 


CdSe 

wz 

48 

50.80 

62 

52 

55 


zb 

41.5 

49.7 

66 

53.4 

45.1®> 


rs 

53 

67.4 


64 


CdTe 

wz 

55 

60 


57.4 



zb 

57 

38.4 

52 

42 

4pl 


rs 

45 

51.86 


53.1 
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ZnO ZnS ZnSeZnTeCdO CdS CdSeCdTe 


FIG. 2. (Color online) Comparison of percentage relative er¬ 
ror in the predictions of (a) lattice constants (b) bulk-moduli 
(c) band-gap (d) energies of the d band calculated using PBE, 
HSE, and ACBNO functionals. All the absolute values are listed 
in Tables The error s are c alculated relative to the ex¬ 

perimental values in Tables mi . In cases where there is more 
than one experimental value, we take the average of the min¬ 
imum and the maximum value as the reference experimental 
value. Comparison is made only for the stable phases of the 
II-VI compounds. Black circles are the PBE values, open red 
squares are the HSE values and open blue diamonds are the 
ACBNO values. The experimental values is referenced by the 
black solid line at 0. 


B. Electronic properties 

The band-structures of the eight semiconductors in all 
the three different crystal phases were computed to assess 
the veracity of the ACBNO functional in the prediction 
of band gaps. For comparison, we computed the band- 
structures within the PBE and the HSE functionals. All 
the band structures were calculated at the theoretical equi¬ 
librium volume optimized within each functional. In the 
wz and zb compounds, the band-gap is direct while in the 
rs compounds the gap is indirect between the L and F 
point. A summary of the values of the energy-gaps {Eg) 
for the different compounds is presented in Table [V| as cal¬ 
culated by the ACBNO, PBE and HSE functionals along 
with references to the SIRC and experimental values when¬ 
ever available. 

Experimentally, the band gap in both the ZnX and the 
CdX series decreases with increasing lattice constants z.e., 
the band gap decreases as we proceed from S to Se to Te. 
Both ZnO and CdO show the anomaly of a lower band-gap 
compared to the ZnS and CdS respectively, in spite of a 
lower lattice constant. This trend of decreasing band gap 
is described by all the functionals: PBE, HSE and ACBNO. 
The absolute value of the band-gap, however, is severely 
underestimated by PBE: in the oxides, especially, we see 
an error bar of more than 60% (See FigJ^(c)). The ACBNO 


functional reduces this error to less than 20%. The band- 
gaps predicted by ACBNO for the eight semiconductors 
studied here lie within 0.5 eV of the experimental value. 

Before discussing the similarities of the HSE and the 
ACBNO functionals, we first look at the general features 
of the calculated band-structures. Representative band- 
structures for the eight semiconductors in their stable 
phases calculated with the ACBNO and the PBE func¬ 
tional are shown in Figures and (band structure of all 
the other phases are compiled in the Supplementary In¬ 
formation). Within PBE, (black lines), the overall band 
profiles are very similar in both the ZnX and CdX series. 
The lowest energy manifolds (below -10 eV) are rather 
narrow and entirely of s character. At higher energy we 
found the bands derived from the d orbitals of Zn and Cd 
that are just below (or little hybridized) with the oxygen 
and chalcogen p states forming the manifolds at the top 
of the valence band. The conduction band is dominated 
by the cation s states. As we proceed from S to Te, the 
chalcogen p states shift up in energy and this results in the 
disentangling of the d bands. Compared to ZnS, in CdS 
the Cd d bands are already well-separated from the S p 
bands. As seen from the PBE band-structures, in both 
ZnO and CdO the d states are more dispersed compared 
to the chalcogenides. The presence of a spurious p-d hy¬ 
bridization in ZnO is a macroscopic manifestation of the 
inadequacy of traditional (semi)local functionals to pre¬ 
dict correctly the electronic properties of materials where 
strong electron localization is present and results in the 
strong underestimation of the band-gap: 0.85 eV versus 
an experimental value of 3.4 eV. Finally, CdO is predicted 
to be a semi-metal by PBE. As it is clearly evident from 
the figures, the ACBNO bands (red lines) maintain their 
overall character as discussed above but display a much 
wider separation between valence and conduction bands. 
The relative position of the d bands is shifted to lower 
energies, thus promoting the disentanglement of the Zn d 
and O p manifolds. For instance, the ACBNO converged 
values for ZnO and CdO energy gaps are 2.91 eV and 0.70 
eV, respectively, in better agreement with the experimen¬ 
tal values of 3.4 eV and 0.9 eV. At the same time, the Zn 
d bands of ZnO are shifted by approximately 4 eV with 
respect to the O p bands, removing the spurious hybridiza¬ 
tion present in PBE. 

Note that the ACBNO follows the HSE functional in pre¬ 
dicting the band-gaps. Both the ACBNO and the hybrid 
functionals reduce the static correlation energy associated 
with the localized orbitals and hence provide a similar 
correction. Indeed, Ivady et.al (Ref. [63|) have recently 
shown that the hybrid exchange correlation potential can 
be rewritten mathematically as an on-site Hubbard po¬ 
tential for a system of localized orbitals, thus providing a 
formal justification for our observations. 

In Fig. id) we have plotted the error bar in the po¬ 
sition of the d manifold of Zn and Cd in PBE, HSE and 
ACBNO. We clearly observe the limits of ACBNO and 
HSE in predicting the energy of the d bands relative to 
the valence band maximum. In PBE and HSE, we see 
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TABLE V. Comparison of the energy-gaps (Eg) in eV for the different compounds using three different functionals. All the PBE 
and HSE values are calculated in this work. In the rs phase, two values are listed. The first one is the direct band gap (T-T) 
while the second one in parenthesis is the indirect band gap between L and T point. Eor the zb and wz phases, the gap is direct 
from T-T. The experiment values are reported whenever available. 


System 

Phase 

PBE 

HSE 

SIRC 

ACBNO 

Experiment 

ZnO 

wz 

0.85 

2.90 

3.4 

2.91 

3.2®, 3.#3 


zb 

0.59 

2.70 


2.74 

3.2';® 


rs 

2.17 (0.89) 

4.35(2.9^ 


4.17(3.09) 

4.5(2.45)P> 

ZnS 

wz 

2.10 

3.42, 3.3^^ 

3.6 

3.31 

3.8d^, 3.91^^ 


zb 

2.23 

3.4^ 

3.6 

3.42 

3.";®, 3.6d^ 


rs 

2.1(0.83) 

3.7(1.31) 


3.52(1.28) 


ZnSe 

wz 

1.27 

2.46 


1.90 

2.87 


zb 

1.18 

2.3d3S,2.4^3 

2.1 

2.041 

2.7(P1 


rs 

1.4(<0) 

2.60(<0) 


2.02(<0) 


ZnTe 

wz 

1.17 

2.22 


1.82 



zb 

1.043 

2.3d^,2.ld^ 

1.4 

1.79 

2.3^ 


rs 

0.6 

1.4 


1.2 


CdO 

wz 

<0 

0.75, 1.1#^ 


1.23 

0.9ll^ 


zb 

< 0 

1.04 


1.35 

2.18-2.23(0.9-1. 08)N^ 


rs 

0.85(< 0) 

2.01(0.8^ 


1.98(0.70) 

CdS 

wz 

1.36 

2.15, 2.0^ 

2.5 

2.4 

2.d^ 


zb 

1.11 

2.13 

2.4 

2.3 

2.5p^_j2.#> 


rs 

1.71(<0) 

2.9(1.4) 


2.95(1.27) 

1.® 

CdSe 

wz 

0.73 

1.77 

1.3 

1.61 

TP 


zb 

0.52 

1.39 

1.4 

1.38 

1.9(P,1.8P 


rs 

<0 

0.69(<0) 


1.4(0.3) 


CdTe 

wz 

0.64 

1.6 


1.33 

1.8 


zb 

0.58 

1.51, 1.5#21,1.6(*S2l 

0.8 

1.43 

1.9^S1,i.61>S21 


rs 

(<0) 

1.4(<0) 


1.3(<0) 



TABLE VI. Average band energies (Ecz) in eV for the eight 
compounds in their stable phases as calculated by PBE, 
ACBNO and HSE functionals. The experimental values are 
taken from Ref. (8] Eor CdO, the value is taken from Ref. 09] 


Compound PBE HSE ACBNO Experiment 


ZnO 

-5.5 

-6.2 

-9.3 

-7.8 

ZnS 

-6.0 

-7.5 

-11.7 

-9.0 

ZnSe 

-6.5 

-7.5 

-12.9 

-9.4 

ZnTe 

-7.0 

-8.0 

-13.8 

-9.8 

CdO 

-5.5 

-6.8 

-9.5 

-8.7 

CdS 

-7.0 

-7.4 

-12.0 

-9.6 

CdSe 

-7.8 

-8.0 

-13.0 

-10.0 

CdTe 

-8.3 

-8.4 

-14.0 

-10.5 


a negative error of around 40% which indicates that the 
position of the d bands are higher in energy compared to 
the experimental values while in ACBNO, the bands are 
pushed much lower, closer to the experimental position. 
The reason for the overestimation of the d band energy 
in HSE calculations is due to the incomplete treatment of 
correlation energy of the localized d electrons. In ACBNO, 
downward shifting of the Sd or 4d bands increase mono- 
tonically with increasing values of t/eff- The rigid shift of 
the bands arises from a singularity due to the filled char¬ 
acter of the d^^ bands and it is implicit in the definition of 
the Hubbard correction. An extensive discussion for the 
case of ZnO is provided in Ref. [30| 


C. Lattice dynamical properties 

In the wide band gap semiconductors, the incorrect de¬ 
scription of the covalency of the cation-anion bond af¬ 
fects the electrostatic properties of the system and, con¬ 
sequently, the phonon distribution and the coupling with 
the external fields. This problem is more pronounced in 
the oxides compared to the chalcogenides. In ZnO the 
underestimated band-gap within PBE functionals leads to 
an overestimate of the high-frequency dielectric constant 
{s^) and, in turn, a large discrepancy in the LO-TO split¬ 
ting. The PBE value of 5.24 is higher compared to the 
experimental value of 3.14 (Ref. l38j) . Similarly, in CdO 
the PBE dielectric constant is 7.13 compared to the exper¬ 
imental value of 5.3. This trend of overestimation of the 
high-frequency dielectric constant is observed in the other 
chalcogenides as well. 

The ACBNO functional provides a proper description of 
the electronic and structural features of the II-VI semi¬ 
conductors and improves the dielectric and vibrational 
properties leading, in most cases, to smaller deviations 
with respect to the experimental data. The calculated 
high frequency dielectric constants (^oo), the Born effec¬ 
tive charges (Z*) and the zone-center phonon frequencies 
are summarized in Tables Ivml andlvnl We did not find 
any references to vibrational properties calculated within 
the HSE method which is computationally very expensive 
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FIG. 3. (Color online) Band structure of ZnX (X=0,S,Se,Te) within the ACBNO (Red lines) method for the stable wz and 
zh phase. The PBE (Black lines) band structure is also shown for comparison. Brillouin zone integration follows the AFLOW 
standard as discussed in Ref. Ell 



ZnO ZnS ZnSe ZnTe 


FIG. 4. (Color online) Band structure of CdX (A=0,S,Se,Te) within the ACBNO (Red lines) method for the stable rs and 
zh phase. The PBE (black lines) band structure is also shown for comparison. Brillouin zone integration follows the AELOW 
standard as discussed in Ref. EH 




CdO 


CdS 




X \VJ^-f- 1 0 \V—L 

CdTe 


especially for calculating the response functions. 

In this section we also show the accuracy of the phonon 
spectrum for six representative compounds: ZnS, ZnSe, 
ZnTe, CdO, CdS, and CdTe. ACBNO phonons frequen¬ 
cies are compared with the PBE and experimental values. 
The dielectric properties and the vibrational spectrum of 
the II-VI semiconductors were calculated using a coupled 
finite-fields/finite differences approach as discussed exten¬ 
sively in Ref. [38l The PBE and ACBNO phonon disper¬ 
sions are reported in Figures and In the Zn chalco- 
genides the acoustic manifold is equally well reproduced 
within the PBE and the ACBNO formalisms (deviations 
of few cm“^ with respect to the available experimental 
data). For the optic manifold, however, ACBNO improves 
the agreement with experiment. This result may be ratio¬ 
nalized in terms of the more accurate description of the 
bonding. In both CdO and CdTe, the ACBNO values are 
in excellent agreement with the experiments and improve 
the PBE description especially for the optic manifolds. In 
CdS both ACBNO and PBE are in good agreement with 
the experiments. 

IV. CONCLUSIONS 

The procedure to compute from first principles the Hub¬ 
bard correction as implemented in the ACBNO functional 


largely improves the electronic structure with respect to 
PBE. The energy band gap in Cd and Zn oxides and 
chalcogenides computed within the ACBNO formalism is 
very close to the value computed with hybrid functional 
at a computational cost of roughly 10%. The ACBNO 
functional generally improves the structural and the vibra¬ 
tional parameters reaching improved agreement with the 
experiments. Both hybrid functionals and ACBNO have 
limited predictive value when considering the the position 
of the occupied d manifolds, at least in ZnS, ZnSe, ZnTe, 
ZnO, CdS, CdSe, CdTe, and CdO that we investigated. 
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FIG. 5. (Color online) Phonon dispersion of ZnX (X=S,Se,Te) within the ACBNO method for the zb phase (red). The PBE 
phonon dispersions (black) is also shown for comparison. Open blue diamonds represent inelastic scattering data from Ref. 1651 
and |66| We show only the directions in the Brillouin zone that have been measured, see the Supplemental Information for the 
full phonon dispersions. 



FIG. 6. (Color online) Calculated phonon dispersion of CdX (X=0,S,Te) within the ACBNO (red lines) and PBE (black lines) 
method for the zb phase. Open blue diamonds represent experiment inelastic neutron scattering data are from Ref. l67ll681 We 
show only the directions in the Brillouin zone that have been measured, see the Supplemental Information for the full phonon 
dispersions. 
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TABLE VIL Calculated frequencies of optic and acoustic 
phonons at specific high-symmetry points in the Brillouin zone. 
We report result computed with ACBNO and PBE functionals 
for all the II-VI semiconductors.Comparison with experiments 
is also reported wherever available. Experimental data are from 
Refs. l68f CdTe) iGTT CdO) l65jZnSe,ZnS) 

System Mode PBE ACBNO Experiment 


ZnS 

TO(r) 

259 

266 

273 


LO(r) 

329 

364 

344 


TO(X) 

293 

301 

312 


LO(X) 

309 

327 

326 


TO(L) 

273 

280 

286 


LO(L) 

319 

337 

333 

ZnSe 

TO(r) 

199 

198 

213 


LO(r) 

240 

260 

254 


TO(X) 

204 

203 

216 


LO(X) 

204 

203 

225 


TO(L) 

201 

200 

222 


LO(L) 

205 

210 

236 

ZnTe 

TO(r) 

170 

172 

175 


LO(r) 

198 

213 

204 


TO(X) 

164 

168 

183 


LO(X) 

174 

182 

181 


TO(L) 

167 

170 

171 


LO(L) 

168 

177 

177 

CdO 

TO(r) 

234 

251 

266 


LO(r) 

410 

462 

465 


TO(X) 

290 

325 

331 

CdTe 

TO(r) 

132 

134 

138 


LO(r) 

158 

171 

167 


TO(X) 

137 

140 



LO(X) 

137 

140 

147 


TO(L) 

133 

137 



LO(L) 

134 

141 

142 


TABLE VIII. High frequency (soo) dielectric constant and Born 
effective charges (Z) of the II-VI compounds computed with 
ACBNO and PBE. The experimental data are from Ref. [65L 

M _ 

z 


System PBE ACBNO Experiment PBE ACBNO Experiment 


ZnS 

5.98 

4.12 

5.13 

2.02 

2.10 

2.1 

ZnSe 

7.40 

4.98 

5.70 

2.12 

2.17 

2.21 

ZnTe 

9.38 

6.01 

7.2 

2.22 

2.19 

2.0 

CdO 

7.13 

3.52 

5.3 

2.40 

1.97 


CdTe 

9.13 

5.56 

7.1 

2.45 

2.34 

2.06 
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